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Abstract 

Dissipation may cause two initially entangled qubits to evolve into a separable state in a finite 
time. This behavior is called entanglement sudden death (ESD). We study to what extent quantum 
error correction can combat ESD. We find that in some cases quantum error correction can delay 
entanglement sudden death but in other cases quantum error correction may cause ESD for states 
that otherwise do not suffer from it. Our analysis also shows that fidelity may not be the best 
measure to compare the efficiency of different error correction codes since the fidelity is not directly 
coupled to a state's remaining entanglement. 

PACS numbers: 03.67.Mn, 03.65.Yz, 03.67.Pp 
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I. INTRODUCTION 



The view of entanglement has gone from being a consequence of quantum theory that 
"no reasonable definition of reality could be expected to permit" [l] to a physical resource 
on which quantum technologies are based Q]. For that reason it is important not only to 
understand entanglement dynamics, but it is desirable to be able to manipulate entanglement 
states in a practical way. 

Entanglement dynamics has been studied in several contexts. When interacting with a 
reservoir, one would naively expect the entanglement between two systems to vanish asymp- 
totically. However, inspired by the entanglement breaking quantum channel theorems [3|, 
Diosi showed that if two systems are initially entangled, and they are subjected indepen- 
dently (or only one of them) to a depolarizing channel, the systems will become disentangled 
in a finite time 4]. This result was subsequently extended for continuous variables in j|. 
Around the same time, Yu and Eberly |6] found that two initially entangled qubits placed 
in their own cavities can become disentangled via spontaneous emission at a rate larser 
than the decoherence rate. For certain initial states, the time of disentanglement is finite, 
a phenomenon termed "entanglement sudden death" (ESD) [7j. Lately, it has been found 
jhat ESD is present in different scenarios for certain classes of states, (see, for example, 



111 ] and the references therein). Moreover, this phenomenon was recently demon- 



strated in an experiment 121 ]. 

A worry, expressed in the literature, is that when an entangled state has suffered from 
entanglement sudden death due, e.g., to dissipation, and become disentangled, no (local) 



recovery operation can restore the state 



ll| . This means that the initial quantum resource 



carried by the state is irrevocably lost. ESD has hence been seen as a potential hurdle, or 
even impasse, for the development of, e.g., large scale quantum computing. 

To combat the decoherence of quantum states, several schemes have been proposed during 
the last decade. When the considered subsystems are coupled to a collective environment, 
one can make use of decoherence- free subspaces 131 ] . Another technique, that does not 
necessarily require that the subsystems are coupled to the same environment, is the open- 
loop control 14J. However, subsequently it was shown that the latter technique is the 



dynamical restoration equivalence of the decoherence-free subspaces 



15 



16]. 



Another option to deal with noisy or dissipative channels is quantum error correction 
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(QEC). In this scheme one can restore a "distorted" state by the introduction of redundancy 



by coding 



3,13, 



19 



2Q|, 
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261 ] . In general, QEC works in three stages: 



At the first stage, every physical qubit is encoded onto a logical qubit (code) composed of 
several physical qubits. The logical qubits then interact with a non-ideal (noisy) channel, 
that "distorts" the state. This distortion is unwanted, but will occur in most realistic 
channels. At the second stage, one performs a measurement to find the error syndrome, 
i.e., what kind of error the channel caused and in which physical qubit it occurred. In the 
third and final stage one restores the distorted logical state into the initial logical state by 
applying a recovery operation, conditioned on the informatio n gi ven by the error syndrome. 
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20| | . and these enables one to 



Several well known error correcting codes exist [13, 
correct a "general" error (e.g., represented by the Pauli operators) in a single qubit. (If two 
or more qubits are erroneous, more sophisticated codes must be used to restore the state.) 
Typically the loss of fidelity after correction is of the order of 7 2 , where 7 is the probability 
of having one error in one of the physical qubits. This result should be compared with the 
fidelity without coding and error correction, where the loss of fidelity scales as 7. 

In open quantum systems not all the possible errors can be represented by the discrete 



Pauli operators ft 



codes like 



17 



lat represent bitflips, phaseflips and both), e.g., amplitude damping. The 
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20j ] are not designed to correct open-system errors. Such errors 



do not fulfil the so-called QEC criteria 



20, 



211 ]. and they can therefore not be exactly 



corrected. To address this problem, a channel-adapted quantum error correction scheme 
has been proposed in 3], where a code is presented that is able to correct a single general 
(discrete Pauli) error and is at the same time perfectly adapted to correct errors due to 
spontaneous emission. Concurrently, it was shown that small deviations from the QEC 
criteria will not change the fidelity qualitatively (in terms of the order of 7), and sufficient 
conditions for an approximate QEC were given [24]. Recently, numerical algorithms to 



3, 3^ 



generate channel-adapted recovery operations using well known codes such as [17J, ll8|, ll9j, |20J 
were presented J25]. In [26] channel-adapted codes for an amplitude damping channel are 
developed, where the encoding is described in the stabilizer formalism. Different recovery 
operations are compared, specially optimal recovery and channel-adapted stabilizer recovery. 
The conclusion from 26] is that the improvement in fidelity by the optimal recovery is small 



compared to the cost of implementation. 

In this paper, we will investigate to what extent ESD can be combatted using quantum 
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coding adapted for dissipative channels. The hope is that the coding, which spreads the 
initial two-qubit entanglement over many qubits, may prevent ESD altogether. Somewhat 
surprisingly, the result is that while coding and error correction will indeed help for small 
amounts of dissipation, quantum coding actually makes matter worse when the dissipation 
is large. In particular, states that do not suffer from ESD when uncoded will, in coded form, 
display ESD. 



II. THE MODEL 



Entanglement sudden death is cased by dissipation. Consequently we will consider an 
amplitude damping channel, that for one qubit can be represented in terms of the Kraus 
operators : 



^l 1 ° , °^ 

v / T Tr 7 / \0 

where 7 is the damping, E is the one-qubit no-jump operator that leaves the ground state 
|0) unchanged, but reduces the probability amplitude of the excited state |1) with the factor 
1 — 7, and Ei represents the jump operator that transforms the state |1) into the state |0) 
with probability 7. 

To investigate to what extent a two qubit state can be "protected" against amplitude 
damping and in particular to ESD, we will use the [4, 1] code introduced in |24|, where the 
notation [4, 1] denotes coding of each initial qubit onto a logical qubit consisting of four 
physical qubits. (We have also tested [5, 1] and [6,2] codes with similar results.) The [4, 1] 
code encodes an arbitrary (single qubit) state \tp) = cos a |1) + sina |0) into the logical state 
\<p) L — cos a \1) L + sin a |0) L , were the codewords |0) L and are given by the four qubit 
code 
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26) 



|o) L = (loooo) + 11111)) /A (i) 

|1) L = (|0011) + 11100)) /V2, (2) 



where here, and below, we will use the notation |0000) = |0) ® |0) ® |0) ® |0). 

For simplicity we will suppose that all the physical qubits have the same damping proba- 
bility 7. We will also assume that the amplitude damping channels are independent, which 
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means that if the amplitude damping channel E acts on the one-qubit density matrix p as 

E(p) = E pEl + EhpEl (3) 

the many-qubit amplitude damping channel will be described by the tensorial products of 
the corresponding Kraus operators, e.g., for two qubits, the channel will be described by 
operators of the form 

E m n = E m (g> E n (4) 

where m,n = 0, 1, and the first(second) operator acts on the first(second) qubit. 

In our assumed channel we have a conditional evolution, because the ground state is 
not affected by the damping whereas the excited state is. Unfortunately, some well-know 
QEC protocols such as those in [17, are not intended for this kind of channels. 



Nevertheless, using the same codewords one can find suitable syndrome measurements and 
recovery operators that do not fulfil the QEC conditions [21J, but the "approximate" ones 



24j. Hence, the [4,1] code in (T1J and ([2]) can only "approximately" correct no jump and 
one qubit jumps, in the sense that the distortion caused by the conditional evolution is of 
the second order in the damping parameter 7 [241 ]. 



As was pointed out in [26|], the states reached from |0) L and |1) L through one-qubit jump- 
ing (i.e., removal of one excitation by the damping) are spanned by orthogonal subspaces. 
Every one of these subspaces are given by two vectors {\Rio) , \Rn)}, i = 0, ... ,4, where 
i = denotes no jump, i = 1, . . . , 4 means a jump of (only) the first physical qubit and so 
on. In total 10 vectors spans the no-jump and one-qubit jump space. For no jump we have 

|i?oo) = |0) L , \R 01 ) = \1) L , (5) 

and for the one-qubit jump it is easy to obtain that 



Rio) = 


10111), 


1^11) = 


|0100), 


R20) = 


|1011), 


\R2l) = 


|1000), 


R30) = 


|1101), 


\R-si) = 


10001), 


-R40) = 


|1110), 


-R41) = 


10010) 



In order to have a full four-qubit basis we can add the vectors: 

\R 50 ) = (10000) - 11111)) /v^, 
5 



R51) = 


(10011) - 11100)) /v^, 


Reo) = 


(|1010) + 10101)) /A 


Rqi) = 


(|1001) + 10110)) /v^, 


R70) = 


(11010) - 10101)) /v^, 


\Rn) = 


(|1001) - |0110)) /v^. 



It is not important how one chooses these additional basis elements since they do not cor- 
respond to any correctable error. 

In this basis, let us introduce the following unitary transformation to make the syndrome 
measurement more transparent: 

7 1 

U = ^2^2\iBm(k))(R ki \, (6) 

fc=0 i=0 

where Bin(/c) is the three-bit, binary representation of k — . . . , 7. In an experimental 
realization, this transformation can be incorporated in the syndrome measurement. 

With the transformation given above one can approximately detect no jump in any physi- 
cal qubit, and one-qubit jumps by measuring the three rightmost qubits in the computational 
basis. If the result is 000 there has been no jump, while the results 001, 010, 011, and 100 
indicate that there have been a jump in qubits 1-4, respectively. Actually, e.g., the result 
010 can also indicate jumps in all of qubits 2, 3, and 4 from the state |0) £ . However, this 
code is not designed to cope with multiple-qubit jumps, so it will act on this syndrome as if 
it was caused by a single qubit jump in qubit 2. 

Once the syndrome is measured we will apply the corresponding recovery operation (sum- 
marized in the table below): If we measure Bin(/c) for k = 0, ...,4, we will preserve 
the possible superposition carried by the first qubit by implementing the transformation 
Rk = Yli=o V) l (iBin(fc)|. If we happen to measure the syndromes Bin(/c) for k — 5, 6, or 
7, we know multiple (but not which) jumps have occurred. Hence we cannot correct the 
state but chose to project the system to the state Il/2, where I L = |0) £ (0| £ + |1) £ (1| £ 
is the identity operator in the codeword space. After correction, the state is always in a 
four-dimensional codeword-space enabling one to compute (or measure) the fidelity to the 
initial state or the concurrence in the standard fashion. 
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TABLE I: The syndromes and corresponding recovery transformations. 



Synd 


rome 


Recovery transformation 


Bin(O) 


= 000 


Rq 


= ELo|i>L<«MO)| 


Bin(l) 


= 001 


Ri 


= lto|i>L<iBin(l)l 


Bin(2) 


= 010 


R 2 


= ZS=o|i>i<*Bin(2)| 


Bin(3) 


= 011 


Rs 


= E-=ol^)L(«in(3)| 


Bin(4) 


= 100 


R4 


= E'=ol*M*Bin(4)| 


Bin(5) 


= 101 




Project onto 


Bin(6) 


= 110 




4/2 


Bin(7) 


= 111 







III. DOES QUANTUM ERROR CORRECTION PREVENT ENTANGLEMENT 
SUDDEN DEATH? 

As pointed out in the introduction, it has been shown that ESD is present in different 



scenarios for some bipartite states 4 



"X-states" 



have been widely studied 



i 0, y, y, 0, 



11 



121 ] . In particularl, the so-called 



9j. They are called X-states because their composite 
density matrix has non-zero elements only in its diagonal and antidiagonal, if expressed in 
the basis {|00) , 1 01) , 1 10) , |11)}. They also have the property that the corresponding density 
matrix preserves the X form when evolving under the action of certain system dynamics, 
and in particular under dissipation. The class contain the Bell-states and the Werner-states 
and therefore merits studying. The class also encompasses both states that succumb to 
ESD and states that does not. To limit our study somewhat, we will work with the Bell-like 
X-states of the form 



\4>o) = cos a 1 11) + e i/3 sina |00) , 
|V>o) = cos a 1 10) + e i/3 sina|01). 

For comparison, we shall also study the separable, non-X-state 



(7) 
(8) 



\ip ) = cos a 1 11) + e l/3 sina 1 10) . 



(9) 
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The state \4>o) will suffer from ESD when | tana| < 1. This was experimentally demonstrated 
in 112 1 



Meanwhile, |^o) does not, and \(fo) certainly does not, as is separable. 
We shall now try to protect the information and entanglement carried by the states (J7|)- 
([§]) from the action of independent amplitude damping channels, and particularly from ESD 
(if present) by encoding each qubit using the codewords ([T]) and ([2]), so that we obtain the 
logical states 

|0) = cosa|ll) L + e i/3 sina|00) L , (10) 
\ip) = cosa|10) L + e i/3 sina|01) L , (11) 
\<p) = cosa|ll) L + e i/3 sina|10) L . (12) 

To characterize to what extent the states are "protected" by coding and error correction we 
will use the fidelity to the initial state and the concurrence. 

After damping and error correction the fidelities between the states f|T0l - f|T2l) and the 
corresponding undamped states, are given (to the second order in 7) by: 

1 2 

= 1-y (7-3cos2a-cos4a) + 0(7 3 ), (13) 

^ = 1 - y (7 - cos 4a) + 0( 7 3 ), (14) 
7 2 

= 1 (17 — 6 cos 2a + cos 4a 

8 

-2cos2/?sin 2 2a) +0(7 3 ). (15) 
Meanwhile, the fidelities between the initial and damped corresponding uncoded states are: 

^-0 O = 1 — 27 cos 2 a + 7 2 cos 2 a, (16) 

•7>o = 1-7: ( 17 ) 

7, 



T vo = 1 — —(11 + 4 cos 2a + cos 4a) 
8 

+-^cos 2 a(3 + 5cos2a) + 0(7 3 ). (18) 
8 



As expected, the loss of fidelity obtained is of the order of 7 2 for the error corrected states 
while it is linear for the uncoded states. The exception is that approaches unity when 
a — > 7r/2. This is to be expected since, in this limit, the state is the trivial, nonevolving, 
ground state. 
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FIG. 1: Plot of the fidelity after recovering for the states \(j>) (dashed), \ip) (dotted), and \<p) 
(continuous) together with the fidelity for the uncoded states |0o) (dash-dotted) and \ipo) (dash- 
double dotted). The parameters are /3 = 0, (a) a = tt/12 and (b) a = 7r/4. (In (b) the curves for 
\<p) and \tp) almost coincide.) 

The expressions for the concurrence C of the states are analytically "messy" , and therefore 
we have chosen to simply plot the exact, numerically results. The plotted fidelities are 
also based on the exact expressions and not on the series expansions f|T3l) - f|T8l) . Since the 
dependence of on f3 is weak, we have, rather arbitrarily, set (3 = when making the plots. 
In Fig. [T]we can compare the fidelities ffT3l) - f|T7j) . In general, the fidelity as a function of the 
damping will be larger for the unentangled state \<p) than for the states \<p) and |?/>). The 
fidelity for the state |0) is larger than the one for \ip) for a < it/A, for a < tt/A the opposite 
happen, being roughly the same when a = 7r/4. 



In Fig. [2] we plot the fidelities and concurrences for the coded state |0) after error 
correction, and the corresponding uncoded state |0o) as a function of the damping. In 
Figure [2] (a) it is seen that the quantum error correction makes the state almost twice as 
tolerant to damping before the entanglement vanishes. However, this is not a general truth, 
and this can be seen in [2] (b) where the coded state suffer from ESD while the uncoded state 
does not (because for the uncoded state C = only as 7 = 1). What is also clear is that for 
small values of 7, quantum error correction increases both the fidelity and the entanglement. 
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FIG. 2: Plot of the fidelity (continuous) and the concurrence (dash-dotted) for the state \(f>) after 
amplitude damping and recovery, together with the fidelity (dashed) and the concurrence (dotted) 
for the uncoded state \(fto). The parameters are (5 = 0, (a) a = tt/12 and (b) a = tt/A. 




3 Q 



-■<. 














\ ; 








X : 








a.z o.s 
Damping y 



(a) 



0.2 0.3 

Damping y 



(b) 



FIG. 3: Plot of the fidelity (continuous) and the concurrence (dash-dotted) for the state \ip) after 
amplitude damping and recovery, together with the fidelity (dashed) and the concurrence (dotted) 
for the uncoded state |V>o)- The parameters are j3 = 0, (a) a = tt/12 and (b) a = 7r/4. 

In Fig. [3] we plot the fidelities and concurrences for the logical state and for the 
corresponding uncoded state |^o). When a = 7r/4, the concurrence and the fidelity for 
the uncoded state \ipo) have the same value, 1 — 7, and the curves in Figure [3](b), are 
not distinguishable. Here again, we see that for large values of 7 the coded state loses 
its entanglement in a finite time, while the uncoded state only looses its entanglement 
asymptotically so that C = only for 7 = 1. However, for small values of 7 both the fidelity 
and the concurrence is increased by the use of quantum error correction, as expected. 
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IV. CONCLUSIONS 



With a few illustrative examples we have studied the possibility of "protecting" two-qubit, 
pure, Bell-like states from the action of amplitude damping. As we have seen, quantum error 
correction succeeds in protecting the state to first order in the dissipation parameter, but 
when the dissipation increases beyond a certain point, quantum error correction often makes 
matters worse and not better. While longer and more complex codes may protect the states 
to higher orders in the dissipation parameter, we conjecture that the change is quantitative 
rather than qualitative compared to simpler codes. In particular, we find that quantum 
error correction will not protect a state from entanglement sudden death. On the contrary, 
states that do not succumb to entanglement sudden death may do so after they are coded. 

Our examples also show what is perhaps obvious, but not much discussed, namely that 
fidelity is insufficient to quantify quantum processes and protocols. In all cases, the fidelity 
decreases asymptotically with increasing imperfections. However, the resource that is often 
of primary interest in quantum information, the entanglement, does not decrease asymp- 
totically. After suffering a finite amount of dissipation the entanglement is gone, while the 
fidelity may still seem comfortably high. The fidelity may thus be deceptive as a figure of 
merit. 

While we have only presented results for the [4, l]-code (applied to the two qubits sep- 



arately), we have also made similar analyses for the [5, l]-code given in [19l . |20j and the 
[6, 2]-code in [26] (that encodes both qubits together). These codes will give marginally bet- 
ter results in the fidelity and concurrence, with similar syndrome measurement and recovery 
operations, but the qualitative behavior, discussed above, remains. However, the [5, l]-code 
is a more complex code and will be more difficult to implement in practice. The [6, 2] code, 
in turn, is a nonlocal code and employs non-local recovery operations. Hence, its implemen- 
tation requires creation of entanglement over the same physical distances as the entangled 
qubits. This is a substantial drawback of the code. 

Our final conclusion is not surprising. In order to combat loss of entanglement due to 
dissipation one must employ quantum error correction. It is important, however, to apply 
the recovery operation while the dissipation is still minor, say, 7 < 0.1. For values of 7 > 0.5 
the entanglement may already have vanished, and will typically have done so if coding for 
quantum error correction has been used. 
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